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Abstract. To study the interaction of star-formation and turbulent molecular cloud structuring, we analyse numer- 
ical models and observations of self-gravitating clouds using the A-variance as statistical measure for structural 
characteristics. In the models we resolve the transition from purely hydrodynamic turbulence to gravitational 
collapse associated with the formation and mass growth of protostellar cores. We compare models of driven and 
freely decaying turbulence with and without magnetic fields. Self-gravitating supersonic turbulence always pro- 
duces a density structure that contains most power on the smallest scales provided by collapsed cores as soon as 
local collapse sets in. This is in contrast to non-self-gravitating hydrodynamic turbulence where the A-variance 
is dominated by large scale structures. 

To detect this effect in star-forming regions observations have to resolve the high density contrast of protostellar 
cores with respect to their ambient molecular cloud. Using the 3mm continuum map of a star-forming cluster 
in Serpens we show that the dust emission traces the full density evolution. On the contrary, the density range 
accessible by molecular line observations is insufficient for this analysis. Only dust emission and dust extinc- 
tion observations are able to to determine the structural parameters of star-forming clouds following the density 
evolution during the gravitational collapse. 

Key words, hydrodynamics - ISM:clouds - ISM:magnetic fields - ISM:kinematics and dynamics - ISM:structure 
- turbulence 



1. Introduction 

Stars form through gravitational collapse in clouds of 
molecular hydrogen. Within these clouds random super- 
sonic gas motions associated with interstellar turbulence 
lead to transient density fluctuations, some of which may 
exceed the threshold of gravitational stability and col- 



lapse to form stars (e.g. 


Elmegreen 1993 


Padoan & Nordlund 1999 


. 



Detailed studies of the hydrodynamic and mag- 
netohydrodynamic collapse in molecular clouds by 



Klessen et al. (2000)| and |Heitsch et al. (2001)| show that 
the location, time scale and efficiency of star formation 
are intimately connected to the stochastic properties of 
the underlying turbulent velocity field. To derive a consis- 
tent picture of the processes from turbulent fragmentation 
to star formation it is important to understand how the 
interplay between supersonic turbulence and self-gravity 
shapes the density and velocity structure of star-forming 
clouds. Thus we need a quantitative measure for the evo- 
lution of the turbulence structure during gravitational col- 
lapse. The density and velocity structure have to be com- 



pared on different scales to follow the process of struc- 
ture transformation from the molecular cloud size to the 
scales of collapsed cores and to understand the influence 
of gravity on the turbulent driving and the overall cloud 
structure. 

A variety of different statistical measures have been 
deployed to analyse the structure of turbulence either 
in observed molecular clouds or in numerical simula- 
tions. Many of them are summarised in the reviews 



by Vazquez-Semadeni (2000) and Ossenkopf et al. (2000) 



One technique especially sensitive to the discrimination of 
the relative structural variation on different spatial scales 
is the A-variance introduced by Stutzki et al. (1998). It 
provides a good separation of noise and observational ar- 
tifacts from the real cloud structure and for isotropic sys- 
tems its slope is directly related to the spectral index of 
the corresponding power spectrum . Bcnsch ct al. (2001) 
and Mac Low fc Ossenkopf (2000) have shown that it can 
be applied in an equivalent way both to observational 
data and hydrodynamic and magneto-hydrodynamic tur- 
bulence simulations allowing a direct comparison. 
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Their investigations, however, neglect the influence of 
gravitational collapse on the structure formation so that 
their conclusions may be limited when applied to star- 
forming regions. It is essential to include the effects of 
self-gravity for the analysis of star-forming regions. It is 
the aim of the current study to close this gap and investi- 
gate the interaction between turbulence and self-gravity. 
We apply the A-variance to characterise the structure in 
numerical models of driven and decaying self-gravitating 
supersonic (magneto-)hydrodynamic turbulence and com- 
pare the results to observed regions of star formation. 

In Sect. | we introduce the molecular cloud models 
analysed here and the numerical methods used for their 
generation. In Sect. [| we apply the A-variance analysis 
and discuss the time evolution of the structure as local 
collapse sets in and protostellar cores form and accrete 
mass. Using various models and comparing the density 
and velocity structure we demonstrate how power is build 
up at the different scales during star formation. Sect. [| 
provides a comparison with observational data including 
dust emission and molecular lines and discusses the result- 
ing differences. We summarise our results in Sect. ||. 



2. Turbulence Models 

The large observed linewidths in molecular clouds im- 
ply the presence of supersonic velocity fields that carry 
enough energy to counterbalance gravity on global scales 
(Williams et al. 2000). As turbulent energy dissipates 
rapidly, i.e. roughly on the free-fall time scale (Mac Low 
et al. 1998, Stone et al. 1998, Padoan & Nordlund 1999), 
either interstellar turbulence is continuously replenished 
in order to prevent or considerably postpone global col- 
lapse, or alternatively, molecular clouds evolve rapidely 
and never reach dynamical equilibrium between kinetic 
energy and self-gravity (Ballesteros-Paredes et al. 1999a, 
Elmegreen 2000). Supersonic turbulence may be driven 
by Galactic rotation and supernovae on large scales and 
by outflows from young stellar objects on small scales (see 
e.g. ic Avillcz 2000] ). The current models investigate both 
cases, continuously replenished and freely decaying turbu- 
lence. 

We select a set of numerical models mostly from pre- 
existing studies that spans a large range of the parame- 
ter space relevant for Galactic molecular clouds. We anal- 
yse the time evolution of their density and velocity struc- 
ture as gravitational collapse progresses. Altogether, we 
include eight models summar ised in Tab. They differ 
in the scale on which turbulent driving occurs, the persis- 
tence of this driving, the inclusion of magnetic fields, and 
the numerical algorithm used to solve the hydrodynamic 
or magnetohydrodynamic equations. 

In the present analysis we neglect feedback effects from 
the produced young stellar objects (like bipolar outflows, 
stellar winds, or ionising radiation from new-born O or B 
stars). Thus our results will hold only for early stages of 
star-forming systems. 



2.1. The driving process 

To represent hydrodynamic turbulence which carries most 
energy on large scales we consider a freely decaying (model 
Sdl) and a driven case (S01). Turbulence that is driven at 
large wavelength leads to a clustered mode of star forma- 
tion (Klessen et al. 2000) and appears to yield the most 
appropriate description of molecular cloud dynamics (see 
e.g. Ossenkopf & Mac Low 2001). Furthermore, we con- 
sider decaying and driven small-scale hydrodynamic tur- 
bulence (models Sd2 and S02, respectively). 

In all driven models the energy input is adjusted 
to yield the same turbulent Mach number M lms — 
10. The Mach number is calculated from the three- 
dimensional rms velocity dispersion, M ms = "rms/cs = 
y / 2Eki n /M /c s , where c s is the sound speed, i?kin the to- 
tal kinetic energy and M the total mass in the system. 
Turbulent equilibrium is obtained applying a non-local 
driving scheme, that inserts energy in a limited range of 
wavenumbers such that the total kinetic energy contained 
in the system remains constant and compensates the grav- 
itational contraction on global scales (Mac Low 1999). 

For the decaying models this Mach number holds at 
times t < 0, before turbulence is allowed to decay. The 
driving stops when gravity is "turned on" at t = 
when the stage of fully established turbulence is reached. 
We also consider a case of large-scale contraction from 
Gaussian density fluctuations without any turbulent sup- 
port (model G). This represents the most extreme case of 
clustered star formation (Klessen 2001). 

2.2. The numerical implementation 

All models are computed either using SPH (smoothed par- 
ticle hydrodynamics) , a particle-based Lagrangian scheme 
to solve the equations of hydrodynamics, or ZEUS, a grid- 
based code. We directly compare the two codes in the case 
of driven turbulence. Here we accompany the SPH mod- 
els S01 and S02 by the corresponding ZEUS models H01 
and H02. This has proven to be a very effective means 
of distinguishing physical results from numerical artifacts 
(Mac Low et al. 1998), as these two approaches practically 
bracket the real dynamical behaviour of interstellar turbu- 
lence. The density enhancements are spatially somewhat 
too small and rigid in the SPH model and too extended in 
the grid-based compuations (discussed in detail by Klessen 
et al. 2000). 

In SPH the fluid is represented by an ensemble of par- 
ticles, where thermodynamical observables are obtained as 
local averages (Benz 1990, Monaghan 1992). The method 
is very flexible and can resolve high density contrasts by 
increasing the particle concentration where needed. We 
use SPH in combination with the special-purpose hard- 
ware device GRAPE (Sugimoto et al. 1990, Ebisuzaki et 
al. 1993), which allows calculations at supercomputer level 
on a normal workstation, and adopt periodic boundary 
conditions (Klessen 1997) focusing on the evolution of a 
sub-region within a larger, globally stable molecular cloud. 
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Table 1. Properties of the considered turbulence models together with the resulting time scales 



model 


description 


kd 


a 


numerical method 


further reference 6 


T 5% 


rd 
Jt {{ 


SOI 


driven HD turbulence 


1 . . 


. 2 


SPH 


Blh in KHM 


0.6 


28% 


Sdl 


decaying HD turbulence 


1 . . 


. 2 


SPH 




0.6 


60% 


S02 


driven HD turbulence 


7.. 


.8 


SPH 




> 5.5 


0.6% 


Sd2 


decaying HD turbulence 


7.. 


.8 


SPH 




2.0 


0.0% 


G 


Gaussian density 






SPH 


X in KB 


1.3 


9% 


HOI 


driven HD turbulence 


1 . . 


.2 


ZEUS 


Vlh in KHM 


0.6 


24% 


H02 


driven HD turbulence 


7.. 


.8 


ZEUS 


mh in KHM 


4.7 


0.5% 


M01 


driven MHD turbulence 


1 . . 


.2 


ZEUS 


S hli in HMK 


1.2 


6% 



"Wavenumber of the original driving 

6 Model names in the original papers: KB - Klessen & Burkert (2000), KHM - Klessen et al. (2000), HMK - Heitsch et al. 
(2001) 

c Time at which 5% of the total mass is accreted onto cores in internal units where — 1.5. 
d Mass fraction that is accreted onto cores after one global free-fall time. 



The simulations use 2 10 5 SPH particles for all tur- 
bulence models except the collapse of the Gaussian den- 
sity fluctuations where we use 5 10 5 SPH particles. Shock- 
generated density fluctuations with masses exceeding the 
local Jeans limit become unstable and collapse. Once their 
central nuclei exceed a density contrast of 10 5 , they are 
identified as protostellar cores and replaced by "sink" 
particles (Bate et al. 1995). These inherit the combined 
masses, linear and "spin" angular momenta of the re- 
placed particles, and have the ability to accrete further 
SPH particles from their infalling gaseous envelopes. This 
technique allows us to follow the dynamical evolution of 
turbulent molecular cloud regions over many global free- 
fall time scales. The numerical scheme and the scaling 
properties of the model are discussed in detail by Klessen 
& Burkert (2000) and Klessen et al. (2000). 

For the grid-based simulations we use the astrophys- 
ical MHD code ZEUS-3D (Clarke 1994). This is a three- 
dimensional version of the code described by Stone & 
Norman (1992a, b) using second-order advection (van Leer 
1977) and von Neumann artificial viscosity introduced to 
resolve shocks. Moreover, structures close to the grid res- 
olution are subject to numerical dissipation as discussed 
by Mac Low & Ossenkopf (2000). 

The simulations were performed on a three- 
dimensional, uniform, Cartesian grid containing 256 3 cells 
and assuming periodic boundary conditions. All quanti- 
ties are scaled to match the corresponding definitions in 
the SPH simulations (see also Heitsch et al. 2001). 

2.3. Magnetic fields 

To test the influence of magnetic fields we consider the 
driven turbulence model which carries most energy on 
large scales and add magnetic fields to the ZEUS simu- 
lations. By comparing the resulting model M01 with the 
hydrodynamic model H01 we get a direct measure for the 
significance of magnetic fields for the interplay of turbu- 
lence and self-gravity in structure formation. 



The magnetic field in this model is selected to be a 
major factor where the ratio between thermal and mag- 
netic pressure (3 = Pth/Pm&g = 8nclp/B 2 = 0.04. With a 
turbulent Mach number of Ai rnls = 10 we find that the 
turbulent velocity dispersion exceeds the Alfven speed by 
a factor of 1.4 so that the structure is still essentially de- 
termined by supersonic turbulence and only modified by 
the magnetic field. The mass in the cloud still exceeds 
the critical mass for a magnetostatically stable cloud by 
a factor 2 (Heitsch et al. 2001) so that the field should 
not prevent gravitational collapse. Cases of sub-Alfvenic 
non-self-gravitating turbulence where the whole structure 
is dominated by the magnetic field have been discussed by 
Ossenkopf & Mac Low (2001). 

2.4. Scaling 

The models presented here are computed in normalised 
units. Throughout this paper we give all size values rel- 
ative to the total cube size, all density values relative to 
the the average density in the cube, and all velocities rel- 
ative to the sound speed. Model G contains 220 thermal 
Jeans masses, wherease all other models have 120 ther- 
mal Jeans masses^. If scaled to mean densities n(H2) = 
10 5 cm -3 , a value typical for star-forming molecular cloud 
regions and a temperature of 11.4 K (i.e. a sound speed 
c s = 0.2 kms -1 ), the mean Jeans mass is one solar mass, 
(Mj) = 1M Q , and the size of cube G is 0.34 pc whereas 
all other models have a size of 0.29 pc. The global free- 
fall time scale, as defined by T ff = (37T/32G) 1 / 2 (p}~ 1/2 
with (p) being the average density, is about 10 5 yr. In nor- 
malised time units it follows rg = 1.5. The simulations 

1 Throughout this paper the spherical definition of the Jeans 
mass is used, Mj = 4/3 7rpAj, with density p and Jeans length 
Aj = (jiTLT/Gp) 1 ' 2 , where G and TZ are the gravitational and 
the gas constant. The mean Jeans mass (Mj) is determined 
from average density in the system (p). An alternative cubic 
definition, Mj = p(2Aj) 3 , would yield a value roughly twice as 
large. 
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cover a density range from n(H2) ~ 100cm -3 in the low- 
est density regions to n(fi2) ~ 10 9 cm~ 3 where collapsing 
protostellar cores are identified and converted into "sink" 
particles in the SPH code. 

In this density regime gas cools very efficiently and it 
is possible to use an effective polytropic equation-of-state 
in the simulations instead of solving the detailed equa- 
tions of radiative transfer. The effective polytropic index 
is typically close to unity, 7 e g < 1, except for densities 
10 5 cm~ 3 < n(ti2) < 10 7 cm -3 , where somewhat smaller 
values of 7off are expected (Spaans & Silk 2000). For sim- 
plicity, we adopt a value of 7 e ff = 1, i.e. an isothermal 
equation of state for all densities in the simulations. As the 
choice of "f c g influences the density contrast in shock com- 
pressed gas, this idealisation may introduce some small 
modifications to the dynamical behaviour compared to 
real cloud systems (see Scalo et al. 1998 or Ballesteros- 
Paredes et al. 1999b for further discussion). 

2.5. Model results 

Supersonic turbulence, even if strong enough to counter- 
balance gravity on global scales will usually provoke local 
collapse. Turbulence establishes a complex network of in- 
teracting shocks, where converging shock fronts generate 
chimps of high density. This density enhancement can be 
large enough for the fluctuations to become gravitation- 
ally unstable and collapse. Long-wavelength turbulence is 
dominated by large-scale shocks which are very efficient in 
sweeping up molecular cloud material, thus creating mas- 
sive coherent structures. When a coherent region reaches 
the critical density for gravitational collapse its mass typi- 
cally exceeds the local Jeans limit by far. Inside the shock 
compressed region, the velocity dispersion is much smaller 
than in the ambient turbulent flow and the situation is 
similar to localised turbulent decay. Quickly a cluster of 
protostellar cores forms. Therefore, all our models with 
zero support, decaying, or large-scale driven turbulence, 
similarly lead to a clustered mode of star formation. The 
efficiency of turbulent fragmentation is reduced if the driv- 
ing wavelength decreases. When energy is inserted mainly 
on small spatial scales, the network of interacting shocks 
is very tightly knit, and protostellar cores form indepen- 
dently of each other at random locations throughout the 
cloud and at random times with low efficiency correspond- 
ing to an isolated mode of star formation (Klessen et al. 
2000, Fig. 1 in Klessen 2001). The effect of magnetic fields 
strongly depends on their strength. Whereas for the sub- 
critical case, i.e. for fields which can support a gravita- 
tionally unstable region magnetostatically (Mouschovias 
& Spitzer 1976) no collapse occurs, supercritical weak 
fields reduce the core growth rate perceptibly, although 
they cannot prevent local collapse (Heitsch et al. 2001). 

To quantify the above we indicate in Tab. |l| the time 
when local collapse has accumulated 5% of the total mass 
in dense cores and also give the core mass fraction after 
one free-fall time scale. Hence, we can directly compare 



the A-variance and other the structural properties of the 
different models cither at equal times or at equivalent de- 
grees of collapse. 

3. The evolution of structural properties 

3.1. The A-variance analysis 

The A-variance analysis was introduced by Stutzki 
et al. (1998) as an averaged wavelet transform 



( [Ziclinsky fe Stutzki 1999[ ) to measure the amount of 
structure present at different scales in an i?-dimensional 
data set. The value of the A-variance at a certain scale is 
computed by convolving the i?-dimensional structure with 
a normalised spherically symmetric down-up-down func- 
tion of the considered size, and measuring the remain- 
ing variance. For two-dimensional structures, like astro- 
physical maps, the filter function is easily visualized as a 
"French hat" wavelet with a positive centre surrounded 
by a negative ring and equal diameters of the centre and 
the annulus. The analysis can be applied in the same way 
in arbitrary dimensions by extending the filter function 
to higher dimensions retaining its radial symmetry and 
adapting the value in the negative part to preserve nor- 
malization (Appendix B of Stutzki et al. 1998). The re- 
sulting A-variance as a function of the filter size measures 
the amount of structural variation on that scale. 

As the convolution with the filter function corresponds 
to a multiplication in Fourier space, we can relate the A- 
variance to the power spectrum of the data set. If the 
structure is characterised by a power law power spectrum 



P(\k\) cx \k\~P 



(1) 



the slope a of the A-variance as a function of lag (filter 
size) is related to the spectral index (3 of the power spec- 
trum by = a + E for < /3 < E + 4. Due to the smooth 
circular filter function the A-variance measures the spec- 
tral index in a way which is more robust with respect to 
edge and gridding effects than the Fourier transform. The 
A-variance provides a clear spatial separation of different 
effects influencing observed structures like noise or a finite 
observational resolution. 

Mac Low fc Ossenkopf (2000)| have shown that we can 
translate the A-variance of a three-dimensional isotropic 
structure into the A-variance of the maps obtained from 
the projections of this structure by rescaling with a factor 
proportional to the lag and an additional small shift. This 
guarantees the preservation of the power spectral index [3 
in projection (Stutzki et al. 1998). As we want to compare 
the simulations to observational data we will always use 
the two-dimensional representation of the A-variance re- 
sults also when applying the analysis directly to the three- 
dimensional data cubes of the simulations. To compute 
the A-variance for our models the SPH density distribu- 
tion is assigned onto a Cartesian grid with 128 3 cells. The 
ZEUS cubes have been analysed in full resolution at 256 3 
as well as degraded to 128 3 for a one-to-one comparison 
with the SPH models. Higher resolution helps to extend 
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the dynamic scale range limited by the periodic boundary 
conditions at the large scale end and the numerical resolu- 
tion of the simulations on the small scale end, but does not 
change the general behaviour of the resulting A-variance. 

3.2. The evolution of the density structure 

3.2.1. The collapse of a Gaussian density field 

Before we investigate the interplay between supersonic 
turbulence and self-gravity, let us consider a system where 
the density and velocity structure is dominated by gravity 
on all scales and at all times. Model G describes the col- 
lapse of Gaussian density fluctuations with initial power 
spectrum P(k) oc k~ 2 and maximum density contrast 
Sp/p « 50. The system is unstable against gravitational 
collapse on all scales and forms a cluster of protostellar 
cores within about two free-fall time scales (Klessen & 
Burkert 2000, 2001). The velocity structure is coupled to 
the density distribution via Poisson's equation and there 
is no contribution from interstellar turbulence. 




Fig. 1. Time evolution of the strength of density fluctua- 
tions as function of their spatial scale measured with the 
A-variance for model G. The density p is given in units of 
the average density in the cube, the lag As in units of the 
cube size, and the time t in internal time units where the 
free-fall time Tg = 1.5. 



The time evolution of the density structure is illus- 
trated in Figure [|. Initially, the A-variance a\ (n) is more 
or less constant (a = 0) on scales As £ 0.02, in agreement 
with the initial power spectrum P(k) oc k~ 2 . The steepen- 
ing below As ss 0.02 is produced by the finite resolution of 
the SPH simulations resulting in the blurring of structures 
at the smallest scales. 

The first changes of the variance o~\{n) are confined 
to small scales. Initial fluctuations with masses below the 
local Jeans limit will quickly smear out by thermal pres- 
sure as the system evolves from purely Gaussian fluctu- 
ations into a hydrodynamically self-consistent state (see 
Appendix B in Klessen & Burkert 2000). As these fluctu- 



ations are by far more numerous than Jean-unstable con- 
tracting ones, the A-variance cr\(n) begins to decrease on 
small scales. However, as the central regions of massive 
Jeans-unstable fluctuations contract to sufficiently high 
densities, o~\{n) increases again. This mainly affects the 
small scales as local collapse modifies the density struc- 
ture on time scales of the local free-fall time. At t = 0.7th 
the first collapsed core is identified and is soon followed by 
others. Altogether 56 dense protostellar cores build up. As 
time advances larger and larger scales exhibit noticeable 
signs of contraction. After about one global free-fall time 
collapse starts to involve all spatial modes in the system 
and the absolute magnitude of the density fluctuations fi- 
nally grows on all scales. As the small scale structure dom- 
inates the density structure we obtain a negative slope in 
the A-variance spectrum. In the final step of the simu- 
lation roughly 30 % of the mass is accumulated in dense 
cores and the slope is about —1.7 indicating that a small 
but significant contribution of large scale structure is still 
present, because an uncorrelated iV-body system of gravi- 
tationally collapsed points would correspond to a slope of 
-2 equivalent to a flat power spectrum P(k) — const. 

The flattening at As > 0.2 is due to periodicity. The 
system is not allowed to collapse freely, it is held up against 
global collapse by periodic boundaries which strongly af- 
fect the evolution of the large-scale modes. The graphs of 
A-variance are not extended beyond effective lags of about 
0.4 as the largest filter that we use is half the cube size 
and we have to apply an average length reduction factor 
of 7r/4 on projection to two dimensions. 

3.2.2. The interaction between gravitation and 
turbulence 

To study the interplay between supersonic turbulence and 
self-gravity, we consider four models of interstellar turbu- 
lence which probe very disparate regions of the relevant 
parameter space. In models S01 and Sdl most of the tur- 
bulent kinetic energy is carried on large scales, whereas 
models S02 and Sd2 involve mainly small-scale turbulent 
modes. In S01 and S02, turbulence is continuously driven 
such that at any moment the overall kinetic energy com- 
pensates the global gravitational energy. In Sdl and Sd2, 
the turbulent energy is allowed to decay freely. 

Figure |^ shows cr^(n) for all four models as func- 
tion of time. In the initial plots one can clearly see 
the dominance of the driving scale as discussed by 



Mac Low & Ossenkopf (2000), The introduction of a ve- 



locity field with a certain scale induces a pronounced peak 
in the density structure at a somewhat smaller scale. Thus 
the curves at t = show for the large scale-driven models 
a power law A-variance from about a third of the cube 
size down to the numerical dissipation scale whereas in 
the small-scale driven model the driving feature at about 
0.07 dominates the structure. 

Star formation is a joint feature of all considered mod- 
els. Like in the evolution of the Gaussian density field 
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Fig. 2. Time evolution of the size distribution of density variations as measured with the A- variance for the four 
models of (a) S01, (b) S02, (c) Sdl, and (d) Sd2. The different depicted times are indicated at the left side of each 
plot. 



the gravitational collapse first modifies only the smallest 
scales, hardly changing the global behaviour. As soon as 
local collapse occurs and the first dense protostellar cores 
form and grow in mass by accretion, they represent the 
main density fluctuations. Their power is concentrated on 
small scales and the A-variance exhibits a negative slope. 
The structure resulting from the collapse is very similar in 
the various models when we compare evolutionary stages 
with about the same mass fraction collapsed onto cores. 

In the initially large-scale driven models a turbulent 
cascade covering all scales is already present from the be- 
ginning shown by the highest a^(n)-values at long scales. 
Within this cascade the number of small-scale fluctuations 
is small compared to models S02, Sd2, and G, and they 
are typically part of a larger structure so that they are 
only weakly dispersed in time. This leads to a monotonous 
growth of the A-variance on all scales. As larger and 
larger regions become gravitationally unstable the con- 
traction comprises increasingly larger scales, and finally 
<j\(n) shifts "upward" on all scales while maintaining a 
fixed slope. This situation is similar in all models with al- 
low for large-scale collapse, i.e. it is only prevented in the 
small scale-driven model S02. 



A different temporal behaviour is visible in Sd2, the 
decaying turbulence which was originally driven at small 
scales. Like in the collapse of the Gaussian density field 
we find in the first steps of the gravitational evolution a 
relative reduction of small-scale structure. This decrease 
is due to the termination of the initial small-scale driv- 
ing resulting in a quick dissipation of the existing fluctu- 
ations by thermal pressure if they arc not Jeans super- 
critical, analogously to the Gaussian collapse case. In the 
next steps we notice the production of structure on larger 
scales. The smoothing of small-scale turbulence combined 
with the onset of self-gravity leads to global streaming 
motions which produce density structures correlated on 
a larger scale. Large-scale structures had been initially 
suppressed by the non-local turbulent driving mechanism. 
After less than one free-fall time a kind of self-sustaining 
inertial cascade with a A-variance slope a of about 0.5 
is build up like in all other decaying models and in the 
large scale driven model. After these initial adjustments 
the first protostellar cores form and we find the same dy- 
namical behaviour. 

The time scale to reach a comparable collapse state 
and the final structure that we reach in the simulations 
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differs between the models, mainly determined by the 
strength of the turbulent driving. The exponent a of the 
A-variance in the collapsed state is -1.3 in the models 
containing continuous driving, -1.5 in the models where 
the turbulence decays during the gravitational collapse, 
and -1.7 for the pure collapse of the Gaussian density 
field. This is understandable, as decreasing turbulent sup- 
port leads to enhanced collapse forming stars in denser 
clusters. The final deviation from the uncorrelated field of 
protostars which has a — —2 is thus a measurable indica- 
tor of the turbulent processes in the cloud during gravi- 
tational collapse. The local collapse produces "point-like" 
high-density cores with small overall filling factor whereas 
most of the volume is supported by hydrodynamic tur- 
bulence. In the large-scale driven model 30 % of the mass 
and in the decaying models 60 % of the total mass has 
turned into cores within the considered time interval. In 
the model continuously driven at small scales only 3 % of 
the mass is in stable cores despite a considerably longer 
simulation time. The small-scale driven turbulence leads 
to the least efficient star formation in an isolated mode, 
whereas the other cases result in the formation of stel- 
lar aggregates and clusters (see also Klessen et al. 2000, 
Klessen 2001). 

Beside the deceleration of collapse the continuation of 
the initial driving has almost no influence on the result- 
ing density structure as soon as the first stable cores have 
formed. It only maintains a constant level of velocity fluc- 
tuations in the main volume of the cloud which is dom- 
inated by low density gas, compared to a homogeneous 
reduction of these fluctuations in the purely decaying case. 

3.2.3. The influence of the numerical model 

Comparing the results of the particle-based SPH and the 
grid-based ZEUS code we can distinguish between numer- 
ical artifacts and physical results, as these two approaches 
practically bracket the real dynamical behaviour of inter- 
stellar turbulence. In Fig. || we compare the A-variance 
plots of the density structure obtained for the driven cases 
using either SPH or ZEUS, at the beginning of the gravita- 
tional collapse and in a step where the structure is already 
dominated by protostellar cores. 

The scaling behaviour of the density structure does not 
differ between both types of simulations but the absolute 
magnitude of the density fluctuations as seen in the total 
value of the A-variance is somewhat larger for all ZEUS 
models. In the first steps of the large-scale driven models 
both numerical approaches still agree approximately but 
during the evolution the scale dependent density varia- 
tions become about a factor five higher in the ZEUS model 
than in the SPH approach. In the small scale driven mod- 
els we can notice a clear difference already at the begin- 
ning of the simulations. This is consistent with the differ- 
ent effective resolution of the methods. Whereas the SPH 
code can provide a very good spatial resolution around 
the collapsing dense regions the general spatial resolution 
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Fig. 3. Comparison of the A-variance measured from the 
particle-based and the grid-based simulations at the be- 
ginning and at about the same timestep of the gravitation- 
ally collapsed state. Plot (a) shows the large-scale driven 
model, plot (b) the small scale-driven situation. 



obtained with 2 10 5 particles is lower than in the ZEUS 
simulations on a 256 3 grid. Thus, the damping of struc- 
tures at small scales due to the finite resolution of the 
code is slightly stronger in the SPH simulations than in 
the corresponding ZEUS models. One can for instance see 
that there is a virtual reduction of structure below 0.01 
which is approximately the radius of the sink particles in 
the SPH code. 

As the SPH resolution is explicitly density dependent 
it is also reduced on all larger scales in low density re- 
gions. This virtually smears out part of the structure on 
all scales. Consequently, the A-variance shows lower val- 
ues on all scales than in the grid-based approach because 
it is not biased towards high-density regions like SPH and 
the observations. The effect is larger in the small-scale 
driven models as the same number of SPH particles has 
to represent more shocks than in large-scale dominated 
cases further reducing the effective resolution. Moreover, 
the resolution worsens during the collapse evolution as 
SPH particles "vanish" in the sink particles. Thus, the 
ZEUS simulations are preferential due to their higher res- 
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olution if one is interested in the absolute value of the 
A-variance whereas they provide no essential advantage 
for the study of the scaling behaviour. 

3.2.4. Magnetic fields 



an additional contribution to the overall isotropic pres- 
sure so that the collapse is somewhat delayed relative to 
the hydrodynamic case but the general structure does not 
deviate from the hydrodynamic simulations. Thus there is 
no need to discuss the weak-field situation here separately. 



3.3. The velocity structure 



0.1 




MHD to 1.0 

HD t= 1 .0 
MHD t = 
HD t = 



Fig. 4. Comparison of the A-variance for the large-scale 
driven model in the hydrodynamic case or the situa- 
tion with strong magnetic fields, both computed with the 
ZEUS code. 



As discussed by pssenkopf fc Mac Low (2001) mag- 
netic fields hardly change the general scaling behaviour in 
interstellar turbulence but create anisotropies in the ve- 
locity field and therefore aligned density structures. Since 
the A-variance cannot measure anisotropies in the density 
structure we do not expect to detect the influence of the 
magnetic field on the turbulence by the present analysis. 
Fig. H shows the A-variance for the initial step and an col- 
lapsed stage in a large-scale driven hydrodynamic model 
and the equivalent MHD model with a strong magnetic 
field. The initial steps are almost identical but we find 
that during collapse the magnetic field effectively helps 
to transfer structure from larger to smaller scales. Thus 
we confirm the more qualitative conclusion of Heitsch 
et al. (2001) that the magnetic field slightly delays col- 
lapse by transferring part of the turbulent kinetic energy 
to smaller scales. The general slope of the A-variance is 
not changed but we obtain somewhat denser and smaller 
cores and somewhat less large-scale correlation at equiva- 
lent timesteps. 

Computing the A-variance for maps projected either in 
the direction of the initial magnetic field or perpendicular 
to it does not show any significant difference in the den- 
sity scaling behaviour as mainly the shape of the collapsed 
regions is influenced, towards spiral-shaped structures, 
which is not measurable with the isotropic A-variance fil- 
ter. 

We have also tested models with a smaller magnetic 
field where the magnetic pressure in the order of the ther- 
mal pressure or lower. Here, we find that the field acts like 




Fig. 5. Evolution of the A-variance of the z-velocity in 
the collapse of the Gaussian density distribution (model 
G). The velocity v z is given here in units of the thermal 
sound speed c s . 

We can apply the A-variance analysis in the same way 
to the velocity structure in the simulations (Ossenkopf & 
Mac Low 2001). Fig. || shows the evolution of one velocity 
component in the collapse of the Gaussian density fluc- 
tuations. In the first steps where we observe a relative 
reduction of small scale density fluctuations we find a bi- 
modal velocity distribution with either very small or very 
large flows. The surplus of small-scale flows just reflects 
the dissipation of the initial small-scale variations by ther- 
mal pressure. When the first stable cores have formed the 
picture changes towards that of a typical shock-dominated 
medium with a slope a = 2 (Ossenkopf & Mac Low 2001) 
as the result of supersonic accretion onto dense cores along 
the emerging filamentary structure (Klessen & Burkert 
2000, 2001). 

In all driven models we see no significant changes in 
the velocity structure during collapse. This is because the 
A-variance is not focused towards the dense cores where 
collapse motions occur, as these have only a small spatial 
filling factor. Instead, most of the volume is occupied by 
tenuous intercore gas with velocity structure that is de- 
termined by turbulent driving. The A-variance therefore 
exhibits the power-law behavior of shock dominated gas. 

The same hold for decaying turbulence as well. To il- 
lustrate that point, Fig. [6] presents the evolution of the 
velocity structure for model Sdl, where we drive the tur- 
bulence initially at large scales and switch off the driving 
during the gravitational collapse. The changes in the ve- 
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Fig. 6. Evolution of the A- variance of the z-velocity com- 
ponent of the decaying model Sdl. 
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Fig. 7. A-variance of all three velocity components for 
model the large-scale driven MHD model at the initial 
step and after one free-falling time. 



locity structure are only minute. The A-variance follows 
the power law of shock-dominated flows throughout the 
entire evolution. Only the total magnitude of the velocity 
fluctuations decreases slightly during the initial decay of 
turbulence. However, after the onset of collapse, when the 
majority of mass is already accumulated in dense cores, 
the magnitude of a a increases again. The A-variance be- 
comes dominated by the shock structure arising from the 
supersonic accretion flows onto individual cores (similar 
to the late stages of model G). The evolution of the small- 
scale decaying turbulence is not plotted separately, as we 
observe the same behavior. Except during the initial phase 
of turbulent decay where the velocity structure still peaks 
on small scales reflecting the smaller driving wavelength 
used to set up the model (see Fig. ||d for the density struc- 
ture), again after t w 1.5 when the initial turbulence is 
sufficiently decayed away the A-variance arrives at the 
power-law behavior of shocked gas. 

As a point of caution, we note here that the veloc- 
ity structure of all SPH models exhibits a numerical arti- 
fact due to the density weighted selection of the size used 
for the averaging kernel. This produces virtual velocity 
streams with relatively large coherence length at low den- 
sities resulting in a A-variance slope a of about 2.8 sig- 
nificantly above the value of 2 characteristic for a shock- 
dominated medium. As the weighting function is given by 
the density of SPH particles, this effect is only visible in 
the velocity structure and not in the density scaling. The 
corresponding grid-based models do not exhibit a similar 
enhancement of large scales in the velocity structure and 
show a A-variance slope of about 2. 

In Fig. [t] we show the three velocity components of 
the MHD model M01 at the same two timesteps like in 
th e density plot in Fig. In c ontrast to the findings 
of Ossenkopf fc Mac Low (2001) for sub-Alfvenic turbu- 
lence, we see no strong anisotropy of the velocity field. 
The velocity structure along the mean magnetic field (z- 
direction) is very similar to the perpendicular directions 
throughout the dynamical evolution of the system and well 



within the statistical fluctuations expected for large-scale 
turbulence. The velocity structure is still determined by 
supersonic turbulence rather than by the magnetic field 
structure because the turbulent rms velocity dispersion 
exceeds the Alfven speed in this model. The local collapse 
and the formation of a cluster of collapsed cores tends to 
make the influence of the magnetic field on the velocity 
structure even weaker. The magnetic field merely decel- 
erates the gravitational collapse and changes the geome- 
try of the collapsing regions as seen in Fig. ^, but hardly 
changes the global velocity structure in this model situa- 
tion. 

Altogether we find that all models that are allowed 
to evolve freely or are driven at large scales exhibit a 
similar velocity scaling behaviour, characteristic of shock- 
dominated media. This is the effect of the undisturbed 
turbulence evolution and the appearance of accretion 
shocks. Both effects lead to remarkably similar properties 
of the velocity A-variance. Any observed deviation from 
this large-range "Larson" behaviour indicative of shock- 
dominated media will hint the presence of additional phys- 
ical phenomena and could provide constraints on the ini- 
tial conditions and the dynamical state of star-forming 
regions. 

4. Comparison with Observations 

4.1. Dust observations 

To compare the simulations of gravitational collapse with 
observational data of collapsed regions we have selected 
the observation of a star-forming cluster in Serpens by 
Tcsti & Sargent (1998) as this region represents a state of 
star formation that may be similar to the outcome of our 
simulations. 

In Fig. H we sh ow the A-variance for t he 3 mm dust 
continuum map of Tcsti fc Sargent (1998) . It reveals an 



increase of the relative amount of structure from small 
scales towards a peak at 7", an intermediate range which 
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4.2. Molecular line observations 



Fig. 8. A-variance of the dust continuum map in Serpens 
taken by Testi & Sargent (1998). The two vertical lines 
represent the limits of the significant range as indicated 
with the observational data. 



can be fitted by a power-law exponent a — —1.2, and a 
decay with a = — 2 indicating the complete lack of large- 
scale structure at lags above 40". The behaviour at largest 
and smallest scales can be under stood when looking at 



the observational base of the map. [Testi fc Sargent (1998) 



give a resolution for their interferometric observations of 
5.5" x 4.3". Consequently we cannot see any structure be- 
low that size. The fact that our peak falls with 7" some- 
what above the 5.5" resolution limit might indicate that 
the CLEAN beam used in the reduction of the interfero- 
metric data is not exactly Gaussian or slightly wider than 
computed. The whole map is taken with an interferometric 
mosaicing technique (see |Testi fc Sargent 2000 for details) 
without a zero-spacing by complementary single dish ob- 
servations. Thus the map cannot contain any structure on 
scales above the single pointing areas determined approx- 
imately by the size of the primary beam of the OVRO 
antennas of 73". This is in agreement with the lack of 
structure indicated by the A-variance slope of -2 at these 
scales. The two limiting sizes are indicated by vertical lines 
in Fig. ^. Thus we may only discuss the range in between 
disregarding other information that is plotted in the inter- 
ferometric map but that can eventually not be obtained 
from the observations. 

The steepening of the A-variance in the intermediate 
size range from a = —1.2 to a = —2 does qualitatively 
agree with the behaviour observed in most collapse simu- 
lations at small scales but does not match any of them 
quantitatively. For a detailed comparison the dynamic 
scale range covered in the simulations is still insufficient 
due to the periodic boundary conditions constraining the 
large scale behaviour. Hence, we can only conclude that 
the collapse models show the same general structure as 
the dust observations, indicating that they represent a re- 
alistic scenario but we cannot yet discriminate between 
different models using the observational data. 



Bcnsch et al. (2001) provided a detailed A-variance analy- 
sis of the density structure traced by observations in differ- 
ent CO isotopes for several molecular clouds with different 
states of star formation including quiescent clouds like the 
Polaris Flare and clouds with violent star formation like 
Orion A. They found for all molecular clouds a density 
structure approximately characterised by a power law A- 
variance, with an exponent in the range 0.5 < a < 1.3. In 
the best studied cloud one smooth curve connects scales 
larger than 10 pc (where turbulence presumably is driven) 
with the dissipation scale at 0.05 pc (where ambipolar dif- 
fusion processes become important). The positive slopes 
indicate that the density structure seen in the CO iso- 
topes is dominated by large-scale modes. This result is 
consistent with purely supersonic turbulence and appears 
independent of the dynamical state of the molecular cloud 
regions studied, i.e. regardless whether the cloud forms 
stars or not. 

This is somewhat surprising, since we expect that the 
density distribution in star-forming regions is dominated 
by the collapsing protostellar cores on small spatial scales. 
The A-variance spectrum therefore should exhibit a neg- 
ative slope as we demonstrate in Sect. The molecular 
line results are also in obvious contradiction to the Serpens 
dust observations discussed above. 

The explanation for the difference is hidden in the 
radiative transfer problem. A discussion of all major as- 
pects of molecular line transfer in turbulent media is be- 
yond the scope of this study and will be provided in a 
separate paper (Ossenkopf in prep.). Here, it is sufficient 
to concentrate on one effect - saturation at large opti- 
cal depths. Molecular lines like the lower transitions of 
13 CO, frequently used to map the density profile of molec- 
ular clouds, become typically optically thick in the cores 
of clouds at densities in the order of 10 5 cm -3 . The ex- 
act value depends on the transition, the spatial configura- 
tion, temperatures, and the geometry of the radiation field 
but one can always assign a typical density range to the 
transition from the optically thin to the optically thick 
regime. This leads to a saturation of the line intensities 
in dense clumps so that the lines do not trace their inter- 
nal structure but rather see clump surfaces. Moreover, the 
molecules tend to freeze out in dense dark regions (Kramer 
et al. 1999) amplifying the effect that the line brightness 
reflects only part of the column density in dense clumps. 

As we do not want to treat the full radiative trans- 
fer problem here, we give only an estimate for the influ- 
ence of optical depth effects by including a saturation limit 
into our computations. Because the simulations are scale- 
free and the typical saturation density varies for different 
molecules and transitions there is no particular density 
value to be used for this limit so that we have to play 
with different values. 

Figure || illustrates the evolution of the large-scale 
driven model shown in Fig. ||a assuming now that all den- 
sities above a certain threshold are invisible so that they 



V. Ossenkopf, R. S. Klessen, F. Heitsch: The structure of self-gravitating clouds 



11 



0.01 






t=0.0 




_ 

m a) - 




t=0.1 




; a 


t = 0.2 






+ 


t=0.3 






□ 


t = 0.4 


^>'^^ 










t = 0.5 


/ /■' ' * s 








//' k .-'Is 




o 


t=0.7 


/// / . / 












T V 


t=1 .6 








* / / / 








& / / 








* / 




- 3 




0.01 0.1 


1 



As 



0.01 








t = 0.0 


b) '_ 






t = 0.1 




A 


t=0.2 






+ 


t=0.3 






□ 


t-0.4 »J 








t = 0.5 -.'O 






o 


t=o.7 






V 


t=1.6 /V 













10~ 3 0.01 0.1 1 

As 



Fig. 9. Time evolution of the A-variance for the large- 
scale driven model SOI when the density structure is as- 
sumed to saturate at densities of 240.0 (a) and 24.0 (b). 



are equal to the value of the saturation density. In the up- 
per plot we have chosen this limit to be the maximum den- 
sity occurring in the original turbulent density distribution 
before gravitational collapse starts. This is 240 times the 
average density, i.e. a relatively large value compared to 
the dynamic range of molecular line observations. In the 
lower graph the saturation limit is reduced by a factor 10. 

Although only a relatively small fraction of the ma- 
terial appears at densities above the limit the influence 
on the A-variance is dramatic. As the collapsing cores 
produce the relative enhancement of small scale structure 
their virtual removal by the saturation results in an almost 
constant A-variance behaviour during the gravitational 
collapse. Even in the very conservative upper plot where 
we assume that all structures occurring in normal inter- 
stellar turbulence are still optically thin and only the col- 
lapsing cores become optically thick the A-variance stays 
at a positive slope during the entire evolution. In the other 
case where optical depth effects are assumed to be more 
important the A-variance remains at a fixed slope of 0.5 
consistent with the molecular line observations of inter- 
stellar clouds. 



Hence, optical depth effects can easily prevent the de- 
tection of gravitational collapse in molecular line obser- 
vations, since they reproduce the A-variance spectrum of 
a turbulent molecular cloud even if collapse has lead al- 
ready to the formation of protostellar cores. Even in mas- 
sively star-forming clouds most molecules used for map- 
ping trace the diffuse density structure between Jeans- 
unstable collapsing protostellar cores. This gas is still 
dominated by interstellar turbulence. The density con- 
trast in star-forming molecular cloud regions simply ex- 
ceeds the density range traceable by molecular transi- 
tions. Protostellar core densities are so high that 13 CO 
at best traces the outer envelope. Therefore both, star- 
forming and quiescent molecular clouds, exhibit very sim- 
ilar molecular line maps. 

It is essential to resolve large density contrasts mea- 
suring the full density structure to study the influence of 
self-gravity and local collapse in star-forming clouds. This 
can be achieved using dust continuum emission. Indeed, 
the 3mm continuum map of the Serpens cluster by Testi 
& Sargent (1998) shows a density structure that is dom- 
inated by small scales as predicted by our collapse sim- 
ulations. The drawback of dust emission observations is 
the inherent convolution of the density structure with the 
unknown temperature profile. Large dust extinction maps 
could circumvent this problem but require long integra- 
tion times at NIR wavelengths to obtain a dense sampling 
with background stars (Lada et al. 1999). 

5. Summary 

In this paper we investigate the time evolution of the 
structural characteristics of self-gravitating supersonically 
turbulent gaseous systems using the A-variance as statis- 
tical measure. We resolve the transition from purely hy- 
drodynamic turbulence to gravitational collapse, the for- 
mation and mass growth of protostellar cores, and com- 
pare different models of driven and freely decaying self- 
gravitating turbulence with and without magnetic fields. 

Contrary to what is observed for purely hydrodynamic 
turbulence, self-gravitating supersonic turbulence yields a 
density structure that contains most power on the small- 
est scales (i.e. in the collapsed objects) as soon as local 
collapse has set in. This happens in all self-gravitating 
turbulence models regardless of the presence or absence 
of magnetic fields. The A-variance o\(n) exhibits a neg- 
ative slope and peaks at small scales as soon as local col- 
lapse produces dense cores. This is in contrast to the case 
of non-self-gravitating hydrodynamic turbulence where 
o\ (n) has a positive slope and the maximum at the largest 
scales. Our results can therefore be used to differentiate 
between different stages of protostellar collapse in star- 
forming molecular clouds and to determine scaling prop- 
erties of the underlying turbulent velocity field. 

The effect of protostellar collapse, however, is not vis- 
ible in molecular line maps of star-forming clouds, as all 
molecules trace only a limited dynamic range of densi- 
ties. The density contrast in star-forming regions is much 
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larger. 12 CO and 13 CO observations, for example, trace 
only the inter-core gas distribution and at best the outer 
parts of individual protostellar cores. Hence, the density 
structure seen in these molecules is indistinguishable for 
star-forming and non-star-forming regions. 

As resolving high density contrasts is the key for de- 
tecting the effect of star formation in the A-variance, we 
propose observations of dust continuum or of the dust ex- 
tinction instead. These techniques do not have the same 
limitations of the dynamic range and are therefore better 
suited to quantitatively study the full density evolution 
during the star-formation process. This is confirmed by a 
first comparison of our models with the 3mm dust contin- 
uum map taken by Testi & Sargent (1998) in Serpens. 

Acknowledgements. We thank Leonardo Testi and Anneila 
Sargent for kindly providing us the FITS file of their 3mm 
continuum map of the Serpens core. Furthermore we are grate- 
ful to J. Ballesteros-Paredes, P. Bodenheimer, A. Burkert, 
C. Kramer, and M.-M. Mac Low for valuable discussions. 
VO was supported by the Deutsche Forschungsgemeinschaft 
through the grant SFB 494B. RSK acknowledges support 
by a Otto-Hahn-Stipcndium of the Max-Planck-Gesellschaft 
and subsidies from a NASA astrophysics theory program 
supporting the joint Center for Star Formation Studies at 
NASA-Ames Research Center, UC Berkeley, and UC Santa 
Cruz. Computations presented here were performed at the 
GRAPE-Cluster at the Max-Planck-Institut fur Astronomie in 
Heidelberg, at the Sterrewacht Leiden, at the Rechenzentrum 
Garching of the Max-Planck-Gesellschaft, and at the National 
Center for Supercomputing Applications (NCSA). ZEUS 
was used by courtesy of the Laboratory for Computational 
Astrophysics at the NCSA. This research has made use of 
NASA's Astrophysics Data System Abstract Service. 

References 

de Avillez M.A. 2000, MNRAS 315, 479 

Ballesteros-Paredes J., Hartmann L., & Vazquez-Semadeni, E. 

1999a, ApJ 527, 285 
Ballesteros-Paredes J., Vazquez-Semadeni E., Scalo J. 1999b, 

ApJ 515, 286 

Bate M.R., Bonnell I.A., Price N.M. 1995, MNRAS 277, 362 
Bensch F., Stutzki J, Ossenkopf V. 2001, A&A 336, 636 
Benz, W. 1990, in: Buchler J.R. (ed), The Numerical Modeling 

of Nonlinear Stellar Pulsations, Kluwer Acad. Publ., 

Dordrecht, p. 269 
Clarke, D. 1994, National Center for Supercomputing 

Applications Technical Report 
Ebisuzaki T., Makino J., Fukushige T., Taiji M., Sugimoto D., 

Ito T., Okumura S.K. 1993, PAS J 45, 269 
Elmegreen B.G. 1993, ApJ 419, L29 
Elmegreen B.G. 2000, ApJ 530, 277 

Heitsch F., Mac Low M.-M., Klessen R.S. 2001, ApJ 547, 280 
Klessen R.S. 1997, MNRAS 292, 11 
Klessen R.S. 2000, ApJ 535, 869 
Klessen R.S. 2001, ApJ 556, 837 
Klessen R.S., Burkert A. 2000, ApJS 128, 287 
Klessen R.S., Burkert A. 2001, ApJ 549, 386 
Klessen R.S., Heitsch F., Mac Low M.-M. 2000, ApJ 535, 887 
Kramer C, Alves J., Lada C.J., Lada E.A., Sievers A., 
Ungerechts H, Walmsley CM. 1999, A&A 342, 257 



Lada C.J., Alves J., Lada E.A. 1999, ApJ 515, 265 

van Leer, B. 1977, J. Comp. Phys. , 23, 276 

Mac Low M.-M. 1999, ApJ 524, 169 

Mac Low M.-M., Ossenkopf V. 2000, A&A 353, 339 

Mac Low M.-M., Klessen R. S., Burkert A., Smith M. D. 1998, 

Phys. Rev. Lett. 80, 2754 
Monaghan J.J. 1992, ARA&A 30, 543 
Mouschovias, T.C. & Spitzer, L. 1976, ApJ 210, 326 
Ossenkopf V., Mac Low M.-M. 2001, A&A submitted 
Ossenkopf V., Bensch F., Mac Low M.-M., Stutzki J. 1999, 

in Ossenkopf V., Stutzki J., Winnewisser G. (eds.), The 

Physics and Chemistry of the Interstellar Medium, GCA- 

Verlag Herdecke, p. 216 
Ossenkopf V., Bensch F., Stutzki J. 2000, in : Gurzadyan V.G., 

Ruffini R. (eds.), The Chaotic Universe, World Sci. , p. 394 
Padoan P. 1995, MNRAS 277, 337 
Padoan P., Nordlund A. 1999, ApJ 526, 279 
Scalo, J.M., Vazquez-Semadeni E., Chappell D., Passot T. 

1998, ApJ 504, 835 
Spaans M., Silk J. 2000, ApJ 538, 115 
Stone J. M., Norman M. L. 1992a, ApJS 80, 753 
Stone J. M., Norman M. L. 1992b, ApJS 80, 791 
Stone J.M., Ostriker E.C., Gammie C.F. 1998, ApJ 508, L99 
Stutzki J., Bensch F., Heithausen A., Ossenkopf V., Zielinsky 

M. 1998, A&A 336, 697 
Sugimoto, D., Chikada, Y., Makino, J., Ito, T., Ebisuzaki, T., 

Umemura, M. 1990, Nature, 345, 33 
Testi L., Sargent A. 1998, ApJ 508, L91 

Testi L., Sargent A. 2000, in Mangum J., Radford S. (eds.), 

Imaging at Radio through Sub-millimeter wavelengths, 

ASP Conf. Series Vol. 217, p. 283 
Vazquez-Semadeni E. 2000, in: Gurzadyan V.G., Ruffini R. 

(eds.), The Chaotic Universe, World Sci. , p. 379 
Williams J.P., Blitz L., McKee C. 2000, in: Mannings V., Boss, 

A. P., Russell, S.S. (eds.), Protostars and Planets IV, p. 97 
Zielinsky M., Stutzki J. 1999, A&A 347, 633 



